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Abstract. The origin of an attractive van der Waals (vdW) interaction be- 
tween two identical neutral atoms is due to an approximately isotropic polar- 
ization, which has been solved exactly if the distance between the atoms is 
very large compared to atomic size. Therefore, the vdW interaction is always 
assumed to be long range by default. Here, we prove (i) the existence of a 
blue-shifting vdW interaction for both long and short distances, and (ii) the 
original vdW theory does not lead to an attractive interaction between two 
strongly interacting atoms with large non-isotropic polarization. 



1. Introduction 

The first van der Waals (vdW) macroscopic model was introduced by Der- 
jaguin pQ, while Lifshitz developed the comprehensive macroscopic vdW theory 
for solid objects by relating the long-range electromagnetic forces to the vdW in- 
teraction within the Green function and quantum field theory formalism [2J. Sub- 
sequently, exact expressions for spherical objects were derived by Langbein within 
the Lifshitz theory [SJ. However, the vdW forces discussed therein are related to the 
interaction between macroscopic solid bodies. On the other hand, the well-known 
exact solution for the semiclassical two-atom system can be found in Ref . [4 , while 
a more comprehensive quantum mechanical results were obtained by Chang [5] via 
the 1/R multipole expansion. We will revisit this expansion later to point out what 
are the conditions that needed to be imposed so that the 1/R multipole expansion 
is valid (given in the Sections, Further analyses I and II). 

Even in this simple system, the attractive vdW interaction has always been 
a long-range type (for weakly interacting liquids), and there is no known result 
whether attractive vdW interaction can even exist at all for atoms polarized non- 
isotropically. Here, we prove that the attractive vdW interaction can indeed blue- 
shift for large non-isotropic (and asymmetric) polarization, and therefore, can be 
made to be significantly relevant even for short distances (shorter than the inter- 
mediate distances considered in Ref. |S]). 

The motivation to revisit the vdW theory is due to its importance in biological 
and biochemical processes such that the blue-shifted vdW interaction could activate 
biochemical reactions [6]. Here, we will attempt to expose the existence of a stronger 
attractive interaction for strongly interacting atoms or molecules in a liquid phase. 
Such a stronger attraction is also important to investigate chemical reactions due 
to blue shifting hydrogen bonds [7] in the presence of directional polarization [5]. 
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Details on this blue-shifting hydrogen bonds are available here [9l [TQl [TT] . Here, we 
prove that the precise microscopic mechanism for the existence of this blue-shifting 
hydrogen bonds can be generalized to all atoms by means of blue-shifting vdW 
attraction. We also show how does this vdW interaction strength changes with 
different types of atoms. Here, we will prove (i) that there is a blue-shifting vdW 
attraction for short distances between two atoms that are asymmetrically polarized, 
and (ii) the original vdW theory does not lead to an attraction between two strongly 
interacting atoms with large non-isotropic polarization. We will address these two 
issues microscopically and analytically. 



2. Renormalized vdW attractive interaction 



The relevant Hamiltonian for the two-atom system (strictly neutral atoms, but 
with polarizable valence electrons) plotted schematically in Fig. [I] is given by [I] 
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where pi j2 are the electrons momenta, and the ground state energy for the above 
Hamiltonian is (taking R 3> r 12 and after some nontrivial substitutions) [J] 
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where m is the electron mass. Note here that ri and r 2 are indeed vectors, but we 
have aligned the atoms such that their electronic polarization are maximum along a 
straight line where ri — > r\ and r 2 — > r 2 . We now prove why this alignment is valid 
for three-dimensional cases (imagine two polarizable spheres) — suppose ri (from 
sphere 1) is aligned along the x axis, so that |ri| = r\ . The induced polarization 
due to ri causes |r 2 | = r 2 (8,(j)) where 8 < n/4, <j) < 7r /4, 8 is the angle in the xy- 
plane and <f> denotes the angle in the xz-pl&ne. It is easy to verify that iff? > tt/A and 
<fi > 7r/4, then there is an additional induced polarization affecting |r 2 | = r 2 (d,<fi). 
If such a case exists, then we just need to take this additional polarization into 
account (that affects r 2 ), and realign r 2 (from sphere 2) with respect to this new 



r[. In this second case, |r^| = r[ causes |r 2 | 



where 8 < 7r/4 and <j) < 7r/4. 



This means that, 8 and <fi remains the same such that 8 g (0, tt/4] and <j> G (0, 7r/4]. 
For one-dimensional cases, 8 = = 4>, while for a two-dimensional case, one has 
8 € (0,7r/4] and <p = 0. □ 

In Eq. (2.2), h is the Planck constant divided by 2ir, R is the distance between 
two neutral and identical atoms. The subscripts 1 and 2 refer to electron 1 and 2 
bounded within nuclei 1 and 2, respectively. The ordinary Hamiltonian denoted by 
Hq consists of two kinetic energy terms, and two semiclassical harmonic-oscillator 
type electron-ion potential terms, and they are associated to two bounded valence 
electrons within their respective nuclei (see Fig. [l]) . The interaction Hamiltonian, 
Hi captures all the Coulomb interaction between the two atoms depicted in Fig. [T] 
Note that Hi does not consider the Coulomb forces within the atom (between the 
valence electrons and their respective nuclei), which have been taken into account 
byH . 



VAN DER WAALS INTERACTION 



3 



Now, one can renormalize the interaction potential constant, k such that |12j 
(2.3) jfe = fcexp[A£], 

where A = (127reo/e 2 )as, in which, cib is the Bohr radius of atomic hydrogen, e and 
eo are the electron charge and the permittivity of free space, respectively. The above 
renormalization procedure is also related (exactly) to the Shankar renormalization 
technique |13l 1141 115) . Here, £ is known as the energy spacing between any two 
energy levels, and it is also known as the the ionization energy. Now the real 
value for £ for a given quantum matter is unknown, and therefore, we will need to 
make use of an approximation. This approximation can be understood from the 
IET-Schrodinger equation j!2) . 



(2.4) ih 



dt 



h 2 

^— V 2 + Viet 
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*(r,t) = ffiET*(r s t) = (£b±0*(r,t), 



where h = h/2n, h denotes the Planck constant and m is the mass of electron, and 
\1/ (r, t) is the time-dependent many-body wave function. The eigenvalue, Eq ± £ is 
the real (true and unique) energy levels for a given quantum system (for a given 
molecule). Here, Eq is a constant because it represents the energy level spacings at 
zero temperature, and in the absence of any external disturbances. On the other 
hand, £ = Emitter m i s the so called ionization energy. Therefore, all we need to do 

. n j i l t - ^quantum 

now is to nnd a way to obtain £„ lattcr . 

Instead of finding a way to calculate ^"ttor™ > which is in itself an unsolved 
problem (since the invention quantum theory), we make use of the ionization en- 
ergy approximation to approximate Cm^ttcr™- This approximation hinges on this 
proportionality, 

/rj r\ ^quantum ^constituent 

Smatter K Satom 

As stated earlier, £m"ttcr m represents the real energy level spacings of a particular 
quantum system, and it can be regarded as a generalized electronic energy gap. This 
means that, in molecules, Cm"™*!™ * s ^ ne energy level spacing between an occupied 
level (in the highest occupied molecular orbital (HOMO)) and an empty level (in 
the lowest unoccupied molecular orbital (LUMO)). Using this approximation, we 
can now predict the changes to Cmattor" 1 if we know £atom ltuent - The values for 
£atom ltucnt can b e readily obtained from the experimental atomic spectra. The 
ionization energy approximation bec ome s exact if the quantum matter are a toms 
or ions (Catter" 1 = Cfon™ from Ec l- l^sj))- The approximation given in Eq. Q 
has been proven — the indirect proof reads, -ffiET^r, t) — (Eq ± Cniattcr m )^'( r ) *) cx - 
(Eq ± ^atom ltuent )^ r ( r i *)) wnue the second formal proof is direct based on logic and 
the excitation probability of electrons and holes within the ionization energy based 
Fermi-Dirac statistics [T2] . 

From here onwards, any variable found to wear a tilde implies a renormalized 
parameter. Subsequently, E can be written in the form (after making use of the 
series, y/1 + x = l + x/2- x 2 /8 + x 3 /16 - . . .) 



m = I 



(2.6) 
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Figure 1 . Two identical atomic system with discreet energy lev- 
els are sketched where ri,2, M and k denote the coordinate for 
the valence electrons, e~ , the ionic mass and the spring (interac- 
tion potential) constant, respectively. The discreet and quantized 
energy levels strictly imply non-free electrons. Here, R is the dis- 
tance between those two nuclei, and M m where m is the mass 
of a valence electron. We set the coordinate such that r\$ and R 
are all on the x-axis where the charges are arranged in a straight 
horizontal line, in the following order, + — H — . Negative charges 
correspond to the valence electrons from the left and right atoms, 
while the positive charges are from the respective nuclei. In the 
text, we also denote the left atom as atom 1, whereas atom 2 is to 
the right. 



in which, we consider only the first three terms by imposing the condition, k exp[A£] 3> 
(e 2 /27re i? 3 ), in other words (e 2 /2ire R 3 ) n k~ n ~ for n > 4 and n e N cvon where 

Neven is the set of even natural numbers. In addition, cDq = 



k/m. Subsequent 



algebraic rearrangements of Eq. (2.6 1, after subtracting Hujq will lead us to 



(2.7) E(0 - hCo 



^Waals(C) 
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This is the standard vdW interaction energy in its renormalized form (due to the 
renormalizing factor (renormalizer for short), exp[— (3/2)A£]), and note the popular 
vdW factor, 1/R 6 . The original (unrenormalized) vdW formula is given in the curly 
bracket. 

The term, Hloq is the renormalized ground state energy (£?(£)) when H i = 
(in the absence of Coulomb interactions between charges from different atoms or 



molecules). One can observe the correct trend from Eq. (2.7), in which, for a given 



R, when £ shoots to infinity (valence electrons with zero polarization), Vv$aais(0 — > 
0, as it should be. In contrast, Vwaais(£) — > maximum, if £ obtains an allowable 
minimum value. Note here that taking £ — is physically not acceptable as this 
means that the valence electrons are not bound to their respective nuclei, and 
they are free. For the same reason, after replacing A = l/kgT where ks is the 
Boltzmann constant [TH], one cannot take the temperature of a system to infinity 
because physically (or chemically), this does not make any sense either, as there 
can be any number of emergent phases due to quantum phase transitions for any 
finite temperature [T?] . 
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3. Stronger vdW interaction due to large non-isotropic 

polarization 



In any case, the above standard attractive interaction (see Eq. (2.7)) is not 
strong enough because the polarization is either small or, if the polarization is large, 
then it is approximately isotropic (see Fig. [2|A)) because the atoms are weakly 
interacting. In the case of strongly interacting atomic-system, a larger induced 
non-isotropic polarization is obtainable (see Fig. [2^13)). Here, Fig. depicts a 

weakly interacting system with an isotropically polarized atom, while Fig.j2^B) and 
(C) depict two types of strongly interacting liquids with non-isotropic polarization. 
To induce such a non-isotropic polarization, we need a stronger interaction between 



these atoms such that the condition, R r is not violated (to make sure Eq. (2.2) 
is still valid). In other words, we need the non-isotropic polarization to be large 
enough (with small £) to give k exp [A£] ps (e 2 /27reo-R 3 ), not k exp [A£] 3> e 2 /2ireoR 3 . 
After making use of this new condition, we can rewrite Eq. (2.2) to obtain (subject 
to the renormalization procedure developed in Ref. 



(3.1) 



E(£) - ftu) = Vwaals(£) 




— Hujq. 



Eq. (3.1 ) can be rearranged to arrive at 
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by using k — (e 2 /2neoR 3 ) interchangeably. 

Obviously, Vwaais(-R) and Vwaais(C) are negative guaranteed by the inequality 
l/y/2— 1 < 0, hence these interactions are always attractive, while the propor- 
tionality, Vwaais(-R) 1/R 3 ^ 2 implies it is of a stronger-type compared to the 
renormalized standard vdW interaction, Vwaais(£) cx 1/i? 6 (see Fig. [lj). It is worth 
noting that the proportionality Vw aa is oc 1/ R 3 ^ 2 is identical to the vdW interaction 
energies derived for two different geometric conditions [THKHj- In particular, the 
relevant geometric conditions are (a) two infinite-length cylinders with radii r% and 
f2, separated by R, and (b) a cylinder (radius r\ and with infinite length) sepa- 
rated from a two-dimensional sheet (negligible thickness but with infinite length 
and width) at a distance R. Contrary to Eq. (3.2), both (a) and (b) must satisfy 
the condition, R <C ri ^, in which, these classical objects are very close to each 



other [HI. In contrast, Eqs. (3.2) and (3.3) are for two spherical quantum objects 
(see Fig. I]) that satisfy conditions, R ^> ri j2 and &exp[A£] ~ (e 2 /2ne R 3 ) (large 
polarization or weakly bound valence electrons). This large polarization will give 
rise to a non-isotropic polarization in the presence of strong interactions between 
atoms. 



However, the imposed condition R 3> rj.,2 does not in any way, implies Eq. (3.3) 
is inapplicable for any R other than R ^> r\ t 2- For example, what Eq. (3.3) tells 
us is that for a given R such that R ^> n^, if £ is large, then one obtains a large 
attractive interaction. Now, the only way to change the value of £ is to change the 
type of atoms because each atom has a unique set of ionization energies or energy 
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Figure 2. We define here the different types of polarizations dis- 
cussed in the text. The sizes and shapes of the sketches are not 
to scale. The nuclei are at the center, while the lines represent 
the energy-levels, similar to Fig. [I] (A) The atom on the left is 
prior to polarization. After intense laser or at a higher temper- 
ature, the atom is polarized isotropically as we have sketched on 
the right-hand side. The two atoms sketched in (A) can also be 
non-isotropically polarized as sketched in (B), if and only if these 
two atoms interact strongly. Therefore, non-isotropic polarization 
is an induced polarization. This is often the case for strongly inter- 
acting liquids such that there exists an effective electron-electron 
Coulomb repulsion leading to a repulsion between these atoms. 
Apart from that, (C) asymmetric polarization can also occur if 
one of the atom is highly polarizable, while the other atom is the 
least polarizable one, which will give rise to an effective attractive 
interaction. Again, this attraction is only valid for strongly inter- 
acting liquids. However, if these two atoms interact weakly, then 
they are approximately isotropically polarized (regardless whether 
the polarization is large or small) , giving rise to the standard vdW 
attractive interaction for R > r± + r%. Here, R is the internuclear 
distance, while r\ and r-i denote the radii of the two isotropically 
polarized atoms. Here, Fig. [2] contains two-dimensional diagrams 
for clarity, but their analyses apply regardless whether a given sys- 
tem is one-, two- or three-dimensional. 



level spacings. Therefore, when we extrapolate this scenario for any R such that 
R > n +T2, then such an extrapolation does not technically and physically change 
the effect of the renormalizer or £. In other words, large £ still gives rise to large 
attractive vdW interaction, regardless whether the atoms are close to each other or 
far apart. 
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Unfortunately, Eq. (3.2) does not provide any new physical insight other than, 
that it is stronger than the original vdW interaction. Whereas, Eq. (3.3 1 is indeed a 
function of the microscopic variable, £, but physically not yet appropriate because 
Eq. (3.3) does not allow large non-isotropic polarization effect — the repulsive inter- 
action between the non-isotropically polarized electrons from two identical atoms 
separated by R, will give rise to an effective repulsion, not an attractive interaction 
(see Fig. |pB)). As a consequence, Eq. ( |3.3[ ) certainly cannot lead to a stronger 
vdW attraction for large non-isotropic polarization, due to the e-e repulsion be- 
tween non-isotropically polarized atoms. 

Moreover, having understood Eq. (2.7), Eq. (3.3) should look downright dis- 
turbing physically since Vwaais(C) oc %/fcexp [(1/2)A£]: note the missing negative 
sign in the renormalizer. In other words, large vdW interaction energy is obtained 
only if the system is least polarizable, that seems to violate Eq. (2.7). But this 
is a false alarm, because if we look closely, it is neither surprising nor physically 
incorrect. The moment we allow large non-isotropic polarization with a new con- 
dition, fcexp [A£] rs (e 2 /27reo-R 3 ), the attractive interaction must now overcome the 
stronger repulsive electron-electron (e-e) Coulomb interaction (V£~° lonlh ) due to 
non-isotropically polarized electrons. Therefore, the attractive interaction is only 
possible if the situation is reversed (compared to Eq. (2.7)), in such a way that 
increasing £ (decreasing polarization) is favored for a maximum vdW attractive 
interaction (Vwaais(£)), as it should be, to avoid the above Coulomb repulsion due 
to non-isotropic electronic polarization. This means that Vwaals(£) ~~ ^ maximum 
when £ obtains a maximum allowable value, or a minimum polarization. 

Here, the minimum allowable polarization implies the minimum value is large 
enough, but not larger, so as not to activate the repulsive e-e Coulomb inter- 
action between the non-isotropically polarized atoms. In other words, we re- 
quire V£~° lomh -I- Vwaais(£) < 0. As such, we can understand the negative sign 
in the renormalizer of Eq. ( |2.7[ ) exists entirely due to the approximation originat- 
ing from the condition, fcexp[A£] 3> (e 2 /27reo-R 3 ), which allows us to invoke the 
solution in series, as stated earlier. In particular, /cexp[A£] 3> (e 2 /2ireoR 3 ) in- 
directly imposes the least polarizability (large £) condition into the system. On 
the other hand, in Eq. ( |3.3[ ), the least polarizability is directly captured via the 
renormalizer, exp [(1/2) A£]. Thus, even though at first we thought that we have 
successfully incorporated the large non-isotropic polarization effect by replacing 
fcexp[A£] > (e 2 /27re i? 3 ) with /cexp[A£] s» (e 2 /27re i? 3 ), but it turns out that we 
are back to square one, with an additional condition, Vq~^ b + V\Vaais(f) < 0. 

As a consequence, we have actually proven that the attractive vdW interaction 
cannot exist for a system consists of two highly non-isotropically polarizable atoms 
(due to strong interaction), for as long as the e-e Coulomb repulsion between these 
non-isotropically polarized atoms is large. In fact, the existence and importance 
°f ^Coulomb between any two highly non-isotropically polarizable entities has been 
analytically proven and discussed in Refs. [2U1 H2 122 ■ Rightly so, the effect of 
non-isotropic polarization, induced by strongly interacting atoms, which gives rise 
to the e-e repulsive interaction (V^~„ lomb ) has been neglected in the standard (for 
both renormalized and original) vdW formula, given in Eq. (2.7). 

Nevertheless, Eq. (3.2), which allows stronger vdW attractive interaction needs 
to be made compatible with Eq. (3.3) in the presence of large non-isotropic po- 
larization. In order to do so, we need to do something to counter the repulsive 
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effect coming from V£~° lomb . We know Eq. (3.3) requires both atoms to be least 
polarizable such that the nuclei of atom 1 (left-hand side (l.h.s) of Fig. [T]) attracts 
the valence electron of atom 2 (right-hand side (r.h.s) of Fig. [I]), and vice versa, 
thus Eq. (3.3) cannot lead to a stronger vdW attraction. 



However, the strength of this attraction can be greatly enhanced if the polariz- 
abilities of these two atoms are made to be strongly asymmetric (see Fig. [2^0)), 
i.e., these atoms are strongly interacting. Meaning, one atom is non-isotropically 
polarized in the presence of another. In particular, if one of the atoms has large £, 
while the other has small £ (£i.h.s < £r.h.s) then this will guarantee the existence of 
^Cou7omb+^Waais(£) < because Vg~u° mb is also negative by definition. In this case, 
the electron-ion (e-ion) Coulomb attraction (VJS"^'^) i s between the polarized va- 
lence electron of atom 1 (easily polarizable) and the nuclei of atom 2 (difficult to 
be polarized). Therefore, to justify the stronger vdW attraction given in Eq. (3.2 ), 



we require asymmetric polarizabilities between two atoms with strong interactions 
(see Fig.||C)). In this case, an asymmetric-polarization can be properly introduced 
into Eq. ( |3.3[ ) to obtain 



(3.4) 



Vwaals(0 = V ( 



ion 
Coulomb 



2^Waals(£)- 



Equation (3.4) is the one that correctly corresponds to Eq. (3.2) in the presence 
of large asymmetric polarization (see Fig. [2](C)). The second term on the r.h.s of 
Eq. (3.4) represents the effect of two identical and weakly interacting atoms coming 
face-to-face (see Fig.JT]), hence the factor 1/2. Both atoms in Fig. [I] are least polar- 
izable (weakly interacting), and they affect each other equally with approximately 
isotropic polarization. While, V^^^ is defined as the interaction between the 
polarized electron of an easily-polarizable atom and the nuclei of the least polariz- 
able nearest-neighbor atom. Therefore, we do not need to include the factor 1/2 
in VqqJJ^j, due to its definition given above, in which the induced non-isotropic 
polarization is unidirectional, from the polarized valence electron of a strongly po- 
larizable atom, to the nuclei of the least polarizable atom. 

In the absence of renormalization, the new condition reads k ps (e 2 /2tt€qR 3 ), 



and one obtains an exact copy of Eq. (3.2) for Vwaais(-R) and 



(3.5) 



Waals 



(k) 



l 



Here, Eq. (3.5) physically correctly captures the same effect as Eq. (3.3). In partic- 



ular, for a given separation R, large k implies small polarization, which is exactly 
identical to large £ that causes small polarization. Hence, two least polarizable 
atoms (weakly interacting) can give rise to the standard vdW attractive interac- 
tion, for instance, when k — ¥ maximum (from Eq. (3.5)), or £ — > maximum (from 



Eq. (3.3)) then both Vwaais(fc) and Vwaais(C) increase to a maximum value due 
to zilch Coulomb repulsion between the polarized electrons (isotropically polarized 
atoms). As stated earlier, Eq. (3.2) is independent of any microscopic variable 



(that can be related to the atomic structure), and as such it is not useful, except 
it correctly predicts the inverse proportionality between Vw aa is and R?/ 2 . On the 
other hand, the changes in k (from Eq. ( |3.5[ )) is arbitrary, for example, we will 
never know how to implement the changes in k accordingly for different atoms in 



a given system. In contrast, and as proven earlier, only Eqs. (3.3) and (3.4) give 



us the option to accurately extract the microscopic information on two competing 
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interactions, between the vdW attraction and the e-e Coulomb repulsion (due to 
non-isotropically polarized atoms) simultaneously, for strongly interacting atoms 
(via £). 



4. Further analyses I 



Equation (3.2) is shown to vary as 1/i? 3 / 2 , and this stronger vdW attraction 
can only be understood by adding an attractive Coulomb interaction (due to asym- 
metric polarization between two strongly interacting atoms) to the standard vdW 
attraction (Eq. (2.7)). Therefore, due to this added Coulomb attraction, one can- 
not obtain Eq. (3.2) using the 1/R expansion since each individual energy term in 
this expansion represents a specific interaction (Coulomb, dipole, quadrupole, and 
so on) without any additional interaction (that gives rise to non-isotropic polariza- 
tion) added to 1/R expansion |S]. Meaning, each specific interaction is determined 
by a specific m and varies as 1/R" where m = 1,2, which is applicable for 
weakly interacting liquids. Recall here that a weakly interacting system means 
there will be no induced non-isotropic polarization, even though these atoms can 
have large isotropic polarization due to intense laser. To understand this point, we 



derive Eq. (2.7) using the second-order perturbation theory (the first-order energy 
is zero) £Q . The renormalized second-order energy 



(4.1) E^ = 



e 

2ne R 3 



2 ~ |(^ 0) (l)ki|4 0) (l))| 2 |(^ 0) (2)k 2 |4 0) (2))| 2 



where (1) and (2) refer to atom 1 and 2 (see Fig. [I]), while the subscripts, s and t 
denote the ground and excited eigenvalues (E st ) and eigenstates (ip a j), respectively. 
Moreover, 



(4.2) 



Hi 



e 2 rir 2 
" 2ne R 3 ' 



after taking R ^> r\i. We first note 



(4.3) 



Ef\l) - Ef\l) + £(°)(2) - El°\2) = 2 



^fuv - ^fr^o 



where (l/2)feI>o > is the ground-state eigenvalue, while (3/2)HQq > is the first 
excited state eigenvalue for the harmonic oscillator-like two atomic system sketched 
in Fig. [Tj Subsequently, using Eq. (4.3) and 

(4.4) |(^ 0) (l)|r 1 |^)(l))| 2 |(^ 0) (2)|r 2 |^)(2))| 2 = [-V ' 



we can derive Eq. (2.7) from Eq. (4.1) following Griffiths [3] — the renormaliza- 
tion procedure introduced here does not disturb the derivation in any way. As 
predicted earlier, E^ oc 1/R 6 , and is always negative determined entirely by 
the eigenvalue-denominator (ground state eigenvalue — excited state eigenvalue = 
— 2/kDo) in Eq. (4.3). The numerator in Eq. (4.1) is positive-definite. Interestingly, 
the same conclusions (without renormalization) can be obtained from Ref. [5]. For 
example, the coefficient Cq [see Eq. (4.3) in Ref. IS]] represents the vdW attrac- 
tion, and is always negative due to e a (k a ) — e a (k' a ) + £f,(fcb) — tb(k' b ) < because 
ea,b( k a,b,k' a .b) < °i M fc a)l < MMI and l e &(^t,)l < l e &0&)l- Here, a and b denote 
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atom 1 and 2, respectively, k a ^ and k' a b represent both the principal and the angu- 
lar momentum quantum numbers such that k a ^ ^ k' a b , while e a ^ are the respective 
eigenvalues. 

To understand the above inequalities, assume a and b are two hydrogen atoms, 
and therefore, e a (k a ) and €f,(kb) are the ground-state eigenvalues for these atomic 
hydrogen, respectively, while € a {k' a ) and tb{k' b ) are the respective excited-state 
eigenvalues for hydrogen atoms. Since, £ , s round = —13.6 eV = e a {k a ) = €b(kb), 
-©excited = —3.4 eV = e a (k' a ) = £b{k' b ), then you can immediately see why the de- 
nominator [in Eq. (4.6) of Ref. [5]] is always negative (as predicted earlier) for as 
long as the polarization is from a lower energy level (strongly bounded electron) to 
a higher energy level (weakly bounded electron). 

Add to that, Cq oc 1/R 6 originates from the interaction potential operator, 
Vz oc 1/R 3 [see Eqs. (2. 13) and (4.3) in Ref. [5J. In our atomic system, we had 



Hi oc — 1/R 3 [see Eq. (4.2)]. Importantly, regardless whether V3 is negative or 
positive, C 6 is guaranteed to be negative due to the eigenvalue-inequalities given 
above, and note that |(- • • | ± V3 1 • • • )| 2 oc 1/R 6 , which is positive-definite. As 
anticipated, V3 is a specific interaction potential (m = 3), which was identified to 
give the vdW attraction without any additional induced non-isotropic polarization 



added to V3. Therefore, similar to Eq. (2.7), V3 excludes the e-e Coulomb repulsion 



in the presence of large non-isotropic polarization, both in the ground and excited 
states or in the mixed states. To those who are still not convinced, read the following 
analyses. 

5. Further analyses II 

The large e-e repulsion due to large non-isotropic polarization (see Fig. [2^B)) 
has been excluded by invoking certain initial conditions used in Ref. [5] to write 
the potential, 



^ V m 



(5.1) V = , 

m— 3 

The initial conditions were R > r a +rb such that most of the charges are distributed 
within an atomic or a molecular radius. This means that the allowable polarization 
is either small or increase isotropically (see Fig. J2^A)). Isotropically increasing po- 
larization implies that the amount of increased polarization is essentially identical 
in all physical directions, namely, in x, y and z axes. Now, these conditions need to 
be imposed so that Eq. (2.14) in Ref. 51 is valid. Rightly so, Chang has correctly 
invoked the above conditions even for large atomic or molecular polarization such 
that the amount of increased polarization is three-dimensionally isotropic. This in 
turn means that the systems studied there were for weakly interacting atoms or 
molecules. If the atoms or molecules are strongly interacting, then any polarization 
that increased isotropically has got to be invalid. In addition, the condition that 
requires the polarization to increase isotropically cannot hold for R w r a + r b . 

Hence, what we have proven earlier is that isotropic polarization is inapplicable 
for R sa r a + Tb (for small internuclear distances) because we need to take the non- 
isotropic polarization effect into account, which has been included in our formalism 



see the first term in Eq. (3.4)]. Moreover, one should also note that the isotropic 



polarization is not applicable for strongly interacting liquids, even if R > r a + r^. 



Thus, Eq. (3.4) is derived for strongly interacting liquids that allow asymmetric 
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polarization for R > r a + r\, (see Fig. |2jC)). In contrast, isotropic polarization can 
only occur for systems that are weakly interacting and R > r a +rf,. It is physically 
incorrect to enforce isotropic polarization for R r a + Tb or for strongly interacting 
liquids. In particular, for strongly interacting liquids with R > r a + r^, there will 
definitely be a significant amount of induced non-isotropic polarization on one atom 
or molecule that is caused by its neighboring atoms or molecules, which then gives 
rise to this strong asymmetric polarization. This is what we have addressed in this 
paper when we derived Eq. (3.4). 

You may correctly wonder if there is an effective e-e Coulomb repulsion be- 
tween two isotropically polarized atoms or molecules. The effective repulsion will 
only be activated when the electrons in both atoms or molecules are polarized non- 
isotropically in the presence of strong interaction between these atoms or molecules, 
and/or when R w r a + r^. For example, atoms or molecules will be isotropically 
polarized in the gas-phase, if these atoms or molecules are weakly interacting in 
the presence of high temperature or intense laser. This isotropic polarization can 
be large. In this case, Eqs. (2.1 1 and (2.2) are applicable where one just starts 
with a new set of radii for atom 1 and 2, and again obtains Eq. (2.7). The new set 
of radii are larger than the original radii because the electrons are now polarized 
isotropically (see Fig. [2]JA)). In fact, one can go on and improve Eq. (2.7) by incor- 
porating other leading terms due to spin-orbit coupling and quadrupole-quadrupole 
interaction as carried out by Chang [5]. But this additional effects cannot induce 
the e-e Coulomb repulsion, which is due to non-isotropic polarization, which oc- 
curs exclusively when there is a strong interaction between atoms or molecules or 
when the average internuclear distance, R « r a + r& (that eventually will induce 
the strong interaction). 

Now, you should not be carried away into thinking that the standard vdW the- 
ory can handle non-isotropically polarized atoms or molecules just because it has 
incorporated other leading terms, namely, 1 /R 5 (quadrupole-quadrupole) and 1 /R 
(magnetic interaction or due to spin-orbit coupling). These terms do improve the 
potential given in Eq. (5.1 1 by incorporating the above stated physical mechanisms, 
and can be used to improve Eq. (4.2), but they (the other leading terms) do not 
cause the e-e Coulomb repulsion to exist. For example, these other leading terms 
are there to improve the equation for Cq or Eq. (2.7) by incorporating the appro- 
priate correction terms such that each term represents a specific physical effect, 
namely, l/R 3 represents magnetic interaction, whereas 1 /R 5 is due to quadrupole- 
quadrupole interaction. In this case, the Coulomb repulsion is not a correction term, 
but originates from the e-e repulsion due to polarized electrons. This means that, 
the first term in Eq. (3.4) is the leading term that represents the asymmetric polar- 
ization, and is responsible for the stronger vdW attraction in strongly interacting 
liquids. In particular, the above e-e Coulomb repulsion, which is due to the large 
non-isotropic polarization, caused by the strongly interacting atoms or molecules, 
can only be overcome by the first term in Eq. (3.4). As a matter of fact, Eq. (2.14) 
in Ref. [5] has been enforced to make sure that the standard vdW theory presented 
in Ref. [5 is for isotropic polarizations, which does not take large non-isotropic 
polarization effect into account. For example, see the spherical harmonics given on 
page 139 in Ref. [3] and Eq. (2.15) in Ref. [5j to observe the expansion for V m that 
is isotropic by definition, which has been enforced in Ref. |5, . For example, the 
standard vdW theory given in Ref. [5] is valid for weakly interacting system and for 
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R > r a +rb such that the polarizations are (assumed to be) approximately isotropic. 
On the other hand, the vdW theory with asymmetric polarization developed here 
is valid for strongly interacting liquids, and for R> r a + r^. 

In summary, a new van der Waals (vdW) interaction of a stronger type is proven 
to exist by generalizing the original vdW interaction, which is also related to Her- 
mansson blue-shifting hydrogen bonds. The generalization includes three types of 
polarizations between two atoms, namely, isotropic, anisotropic and asymmetric 
polarizations. We show that an isotropic polarization is only valid for R ^> 2r 
(r = T a = rb), which is a well known result within the original vdW theory. But 
the other two polarizations (anisotropic and asymmetric) are shown to be valid for 
all R. Recall here that R is the internuclear distance between two atoms, while r 
is the atomic radius for two identical atoms. Using the above generalization, we 
went on to explain why and how the original vdW attraction does not lead to an 
attractive interaction between two highly polarizable atoms due to large Coulomb 
repulsion between polarized electrons. This explains why i? > r is a necessary 
condition for systems with isotropic polarizations, in accordance with the original 
vdW formalism. However, for R rts 2r, the two atomic system should be considered 
as strongly interacting, giving rise to an anisotropic or an asymmetric polarization. 
In these cases, there is no requirement on R such that these polarizations are valid 
for any R. Hence, these two (anisotropic and asymmetric) polarizations, which are 
not captured by the standard vdW, have been proven to exist, which could lead us 
to understand the (bio)chemical reactions down to an electronic level. 

6. Conclusions 

We have formally proven the limitation in the original vdW formalism, in which, 
the Coulomb repulsion between polarized electrons has been neglected. In particu- 
lar, the Coulomb repulsion is deliberately avoided by requiring small polarization, 
which has been implicitly embedded into the original formalism by requiring the 
interaction potential constant (fc) to be much greater than the Coulomb interaction 
between atoms (molecules). Or, by imposing a large isotropic polarization by in- 
voking spherical harmonics. After taking the large non-isotropic polarization effect 
into account, we managed to obtain a new vdW formula of a stronger type, which 
is proportional to l/i? 3 / 2 . However, the new formalism also leads to a formula 
that is proportional to cxp[A£] where £ is the energy level spacing of a given mole- 
cule or system. Here, large £ means large k, in accordance with the original vdW 
theory. Therefore, large non-isotropic polarization does not favor large attractive 
interaction at all, because of the e-e Coulomb repulsion effect (due to polarized 
electrons) . 

To overcome this intrinsic problem, large asymmetric polarization between two 
interacting molecules was invoked to prove the existence of a blue-shifting vdW 
interaction, in which, the above repulsion has been physically converted into a 
Coulomb attraction via an asymmetric polarization. For example, the attraction 
between the polarized electrons (of an easily non-isotropically polarizable; atom) and 
a nearest neighbor nucleus (of the least polarizable atom) gives rise to this asym- 
metric polarization. This means that, the large vdW interaction (due to 1/i? 3 / 2 ) 
can be justified with the large asymmetric polarization effect. Therefore, there 
exists a blue-shifting attractive vdW interaction between two strongly interacting 
atoms if they are polarized asymmetrically such that one of the atom or molecule 
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is non-isotropically polarized due to interaction with its nearest least-polarizable 
neighbor. 
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